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W. P. ORRICK 

Abstract. We study the question of finding the maximal determinant of matrices of 
odd order with entries ±1. The most general upper bound on the maximal determi- 
nant, due to Barba, can only be achieved when the order is the sum of two consecutive 
squares. It is conjectured that the bound is always attained in such cases. Apart from 
these, only in orders 3, 7, 9, 11, 17 and 21 has the maximal value been established. In 
this paper we confirm the results for these orders, and add order 15 to the list. We follow 
previous authors in exhaustively searching for candidate Gram matrices having deter- 
minant greater than or equal to the square of a known lower bound on the maximum. 
We then attempt to decompose each candidate as the product of a {—1, l}-matrix and 
its transpose. For order 15 we find four candidates, all of Ehlich block form, two having 
determinant (105 • 3''' • 2^**)^ and the others determinant (108 • 3^ • i^^f. One of the former 
decomposes (in an essentially unique way) while the remaining three do not. This result 
proves a conjecture made independently by W. D. Smith and J. H. E. Cohn. We also use 
our method to compute improved upper bounds on the maximal determinant in orders 
29, 33, and 37, and to establish the range of the determinant function of { — 1, l}-matrices 
in orders 9 and 11. 



1. Introduction 

Under the conditions that its elements are real and bounded in magnitude, what is the 
largest determinant that a square matrix can have? By rescaling we may fix the bound 
on the magnitude of the elements to 1, and because of linearity of the determinant in 
its rows and columns, we lose nothing by restricting the entries to the set {—1,1}. For 
Rn G {n X n matrices with entries ±1}, Hadamard [Hij showed that 

deti?„<n"/2 (1.1) 

and that equality can only hold when n = l,2orn = mod 4. It is not known whether 
equality always does hold in such cases, but Paley |Paj conjectured that it does, and 
the lowest order for which the question is unresolved is n = 428. For other orders more 
stringent bounds have been established. Barba |Baj showed that when n is odd the bound 
is 

deti?„<y2^r^(n-l)("~^)/^=5(n). (1.2) 
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This bound is an integer only when n is the sum of two consecutive squares, and it is 
beheved that equahty holds whenever this is the case. However, even for orders as low 
as 85 = 6^ + 7^ this has not been proved. When n = 2 mod 4, the bound, derived 
independently by Ehlich |Ehlj and Wojtas JWo' , is 

deti?„ < 2(n- l)(n-2)("-2)/2. (1.3) 



The bound has been shown to be attainable only if ri — 1 is the sum of two squares |Ehlj . 
and again it is believed to hold in all such cases. The lowest order for which the question 
has not been settled is n = 138 |FKSj . 

We say that a bound is tight when it is achieved infinitely often. All three of the above 
bounds have been shown constructively to be tight \Bt[ |Sp[ fKKSj . The bound (|1.2p is 
actually tight only for n = 1 mod 4. For n = 3 mod 4 Ehlich |Eh2j proved the stricter 
bound 



det Rn<{n- 3)("-^)/^(n - 3 + 4r)"/2(^ + i + Ary/^Ji Hr + I) 

where s = 5 for n = 7, s = 5 or 6 for n = 11, s = 6 for 15 < n < 59 and s=7 for n > 63 
and where r = [^J , n = rs + f and u = s — v. This bound is not integral for n < 59. For 
n > 63 Cohn |Co2j has shown that it is integral only when n = 112t^ ± 28t + 7 for some 
integer t, but even in these cases, it is not known whether the bound is ever attained. 

Denote by md(n) the maximal determinant attained by any member of the set of 
n X n { — 1, l}-matrices and define md(n) := md(n)/2"~^. The latter is the "normalized" 
maximal determinant function where the factor 2""^ common to all nxn {—1, l}-matrices 
has been removed. 

When a determinant is found that attains the relevant one of the above bounds, it is 
immediate that md(n) is just the bound itself, but when the upper bound is not attained, 
finding md(n) can be exceedingly difficult. To date, the only orders in this category for 
which md(n) is known are 3, 7 jWij . 9 |EZj . 11 (Ehlich, as reported in |GKj ). 17 |MKj 



and 21 |C;KMj . This leaves, for n < 30, orders 15, 19, 22, 23, 27, and 29 unresolved. The 
sequence of maximal determinants for 1 < n < 14 is 



n 

md(n) 



1 2 3 4 5 6 7 8 9 10 11 12 13 14 
1 1 1 2 3 5 9 32 56 144 320 1458 3645 9477 



In this paper, we settle the case n = 15. This means the above sequence can be extended 
up to n = 18: 

n 15 16 17 18 

^{n) 25515 131072 327680 1114112 ' 

Our method of proof follows that of Ehlich (reported in |GKj ). Moyssiadis and Kou- 
nias JMKI and Chadjipantelis, Kounias and Moyssiadis |CKMj . and relies on an intensive 
computer search. Compared with all other orders in which md(n) has been determined, a 
considerably greater (but still modest) amount of computer power is required for n = 15. 
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The value 25515 was conjectured to be the maximum by Warren D. Smith |Smj and 
independently by John H. E. Cohn |Coll ICo2] . 

In addition to this result for n = 15, we confirm the known values of md(n) for all 
odd n < 21 excluding, of course, n = 19. We are also able to report improved upper 
bounds on the maximal determinants for matrices of orders 29, 33, and 37. Furthermore, 
we establish the range of the determinant function for matrices of orders 9 and 11. 

2. Some properties of matrices 

While we are here mainly interested in the case n = 15, we also would like to confirm 
existing results in other orders and to say as much as possible about other open cases. 
On the other hand, even orders (the lowest interesting case is n = 22) require a somewhat 
different approach, owing primarily to the absence of a canonical normalization (see be- 
low). Hence we will describe an algorithm that works for odd orders generally. For the 
remainder of this paper we take n to be an odd positive number. 

2.1. Equivalences of { — 1, l}-matrices. Negating some set of rows and some set of 
columns of a matrix changes its determinant by at most a sign. Two matrices related 
to each other by such negations are considered to be equivalent. It is advantageous to 
restrict our attention to a single canonical representative of each equivalence class: 

Definition. A vector of odd length with elements in the set { — 1, 1} is parity normalized 
if it has an even number of positive elements. A matrix of odd order with elements in 
{— 1, 1} is parity normalized all its rows and columns are. 

Lemma 2.1. Any {—1, l}-matrix of odd order may be converted to a unique parity nor- 
malized matrix by a series of negations of rows and columns. 

Proof. Begin by negating all rows containing an even number of entries —1. Since the 
number of rows is odd, the matrix now contains an odd number of negative entries. Thus 
the number of columns containing an even number of negative entries must be even. 
Negating these gives the desired matrix, which is clearly unique. □ 

Permuting a set of rows and a set of columns of a matrix also preserves its determinant 
up to a sign. 

Definition. Two { — 1, l}-matrices R and S are Hadamard equivalent if S* = PRQ for 
some pair of signed permutation matrices (P, Q). 

Note that we have omitted a third determinant preserving operation, transposition, in 
our notion of equivalence. Henceforth, we take { — 1, l}-matrices to be parity normalized. 
On the other hand, transforming a given matrix into a canonical representative of the 
same class under full Hadamard equivalence is computationally expensive, so we do not 
attempt to impose a canonical ordering on the rows and columns of a matrix. 
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2.2. Gram matrices. Let Mr{R) be the Gram matrix of tlie rows of R, 

Mr{R) = RR^. (2.1) 

Mr{R) is unclianged under negation and permutation of the columns of R. Likewise let 
Mc{R) be the Gram matrix of the columns of i?, 

Me(i?) = R^R. (2.2) 

Mc{R) is unchanged under negation or permutation of the rows of R. Both Mr{R) and 
Mc{R) are symmetric, have diagonal elements equal to n, and, assuming R nonsingular, 
are positive definite. They have the same determinant, which must be a perfect square. 
Furthermore, they have identical characteristic equations. This follows from the following 
lemma. 

Lemma 2.2. Let A and B be square matrices. Then det{AB — XI) — det{BA — XI). 

Proof. If B is invertible, the result follows immediately by multiplying the right-hand side 
of the equation on the left by det B and on the right by detB~^, and then combining 

determinants. If B is not invertible, consider B' = B — xl. Since B' is invertible for all 
but finitely many x, we have, in the generic case, that dct{AB' — XI) = dct{B'A — XI). As 
this is a polynomial identity in x, it must be true for all x and for x = in particular. □ 

Since Gram matrices are symmetric, we will often have occasion to talk about operations 
that act on rows and columns simultaneously in the same way. We will therefore refer to 
such operations as, for example, permuting the indices of M. 

When R is parity normalized, the elements of Mr{R) and Mc{R) are necessarily con- 
gruent to n mod 4. Under permutation of the rows of R, the indices of Mr{R) undergo 
the same permutation. Under permutation of the columns of R, the indices of Mc{R) 
undergo the same permutation. 

The search method we employ requires the generation of a set of candidate matrices 
which are potentially Gram matrices of { — 1, l}-matrices. We are interested in { — 1,1}- 
matriccs whose determinant equals or exceeds some threshold, rfmin in magnitude. For a 
given n, this threshold is usually chosen to be the best known lower bound on md(n). To 
be a candidate, a matrix M must possess the properties of Mr{R) and Mc{R) described 
above and must have determinant greater than or equal to d'^^^. The set of such candidates 
is the set A4^(ciniin) whose definition is made precise below. 

Definition. We denote by M.n,p be the set of p x p matrices M having the following 
properties: 

(1) M is symmetric and positive definite; 

(2) the diagonal elements of M are all equal to n; 

(3) the elements of M are integers congruent to n mod 4. 

The set A4n,p{dinin) is the subset of M.n,p whose elements satisfy the additional requirement 

(4) det M>dl,,. 
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We define A4n = ■Mn,n and A1„((imin) = ■A4n,n{dmin) ■ Finally, let A1„ ((imin) be the subset 
of Aln (c^min) whosB elements satisfy 
(5) det M — perfect square. 

Positive definiteness implies that the off-diagonal elements of M are less than n in 
absolute value. 

Just as for { — 1, l}-matriccs, there is a notion of equivalence for candidate Gram ma- 
trices. Negation will not be considered because the congruence property modulo 4 that 
we have imposed on the matrix elements already fixes the signs. 

Definition. Two matrices, Mi,M2 G A4n,p, are equivalent if one can be obtained from 
the other by a permutation of indices. The equivalence class of M will be denoted £m- 

Definition. An element of a matrix M e A4n,p is minimal if its magnitude is the smallest 
allowed value. If n = 1 mod 4 the minimal element is 1; if n = 3 mod 4 the minimal 
element is —1. The minimal vector of length k is the vector of length k, all of whose 
elements arc minimal. 

Definition. A set of indices A C {1. ... .p} is a block of M G Ain,p if, for every pair 
with i & A and j G {1. . . . ,p} \ A, the matrix element M^j is minimal, and if no proper 
subset of A has this property. 

Remark. The blocks of any matrix M G M.n,p are disjoint, and there is a unique decom- 
position of the full set of indices {1, . . . ,p} into blocks. 

Definition. A set of indices A C {1, . . . ,p} is contiguous if for any i,k & A and j G 
{!,..., n}, i < j < k implies j G A. 

Definition. A matrix M G M.n,p has block form if all of its blocks are contiguous. 

We will also use the term "block" to refer to the sub-matrix obtained by removing all 
rows and columns whose indices are not in the block. 

Definition. The partial row Pr,k{M) with M G A4n,p and 1 < r. A; < p is the k- 
dimensional vector formed by taking the first k components of row r of M. 

Notation. For integers a and b we use a <a & to indicate that \a\ < \b\, with similar 
notation for other inequality signs. For integer vectors u and v, u <a v indicates that u 
is lexicographically less than v where comparison of individual vector elements is made 
using <a. 

Definition. A matrix M G A^n,p is lexicographically ordered if its rows are in descending 
lexicographic order, according to the ordering <a. 

Since each row of M is dominated by its diagonal element, this definition is equivalent 
to the criterion that, for all pairs of indices 1 < i < j < p, the partial rows satisfy 

Pi,i-l{M) >aP,-^-l(M). 

Lemma 2.3. A lexicographically ordered matrix M G M.n,p) has block form. 
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Proof. Assume the contrary. Then there is at least one pair of blocks, {A,B), whose 
indices are interleaved. Assume the lowest index in A to be less than the lowest index 
j in B. Then A is partitioned into two subsets: A^j of indices less than j and A^j of 
indices greater than j. There exists i G A^j and k G A^j such that Mj ^ is non-minimal, 
since otherwise A^j and A^j would be separate blocks. 

By the lexicographic ordering of M we have Pjj-i{M) >a Pk,j~i{M). However, all 
entries of Pj j_i(M) are minimal since j is the smallest index in its block, whereas Mk^i = 
Mi^k is non-minimal — a contradiction. □ 

Notation. Let Mi G Mn.p^ and M2 G Mn,p2 We say that Mi <a M2 if either 

(1) there exists an index k such that 1 <k < min(pi,p2)5 = -Pj,j-i(^2) for 
all i < k, and Pfc,fe_i(Mi) <a Pk,k-iiM2), or 

(2) pi < p2 and Pi,i-i{Mi) = Pi^i.iiM^) for all i < p,. 

Clearly the blocks of a lexicographically ordered matrix (considered as sub-matrices) 
are themselves lexicographically ordered. 

There may be many lexicographically ordered matrices equivalent to a given matrix. 
Because <a is a total ordering and A4n,p is finite, any subset of J^n,p has a greatest 
element . 

Definition. A matrix M G Ain,p is lexicographically maximal if it is the greatest element 
of Sm under the ordering <ay. 

Trivially, a matrix that is lexicographically maximal is lexicographically ordered. The 
proof of the following is straightforward. 

Lemma 2.4. A matrix M G Ain.p is lexicographically maximal if and only if all its blocks 
(considered as sub-matrices) are lexicographically maximal and are arranged in descending 
lexicographic order along the diagonal. 

In Section 121 we present a backtracking algorithm to generate a list containing a canon- 
ical representative of each equivalence class in A^nlc^mm)- The canonical representative 
of an equivalence class is taken to be its lexicographically maximal element. The algo- 
rithm builds up matrices by augmenting the leading sub-matrix, one row and column per 
iteration, and is so designed that any sub-matrix M G A4n,p, constructed at stage p, is 
automatically in block form with lexicographically ordered blocks in descending lexico- 
graphic order along the diagonal. To ensure that only canonical matrices are produced 
in the end, there is a check, performed at the completion of any block, that the block 
is lexicographically maximal. In the event that it is not, the sub-matrix is discarded. 
To carry out this check, the following simple procedure, while not optimal, proved to be 
adequate. 

2.3. Procedure IsLexMax. The procedure returns True if the block cannot be trans- 
formed into some lexicographically greater block by permutation of indices, and returns 
False otherwise. It is perhaps simpler to consider how to put a block into its lexico- 
graphically maximal order. There is at least one permutation that accomplishes this. If 
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the set of all such permutations includes the identity, then our routine is to return True, 
otherwise False. 

Let the block be of size p > 2. The procedure goes as follows: we first find the set, V2, of 
permutations that produce the lexicographically maximal 2x2 leading sub-matrix. Since 
this is accomplished by placing a maximal off-diagonal element into the (1,2) position, the 
set V2 can contain permutations of three types, depending on the initial locations of the 
maximal elements: the identity, if a maximal element already was in the (1,2) position; 
a single swap, if a maximal element occurred in a position with an index (row or column) 
1 or 2; a pair of swaps, if a maximal element occurred in a position with neither index 
equal to 1 or 2. 

Now, given the set Vk of permutations that produce the maximal k x k leading sub- 
matrix, 2<A;<p — 1, we produce the set Vk+i- We do this by composing each per- 
mutation in Vk with each single swap {k + 1, j), j > k + 1, that produces the maximal 
{k + 1) X (k + 1) leading sub-matrix (if any exists), and with the identity, if the leading 
[k + 1) X [k + 1) sub-matrix already is maximal. The permutations so produced are the 
elements of Vk+i- It is in fact unnecessary to consider all possible swaps, since some 
indices in the set {k + 2, . . . ,p} may be equivalent (under permutation of indices) to each 
other due to the structure of the matrix. Thus we need only consider one index from each 
equivalence class. 

This process is iterated until either 

(1) a set Pfc does not contain the identity, in which case we terminate the procedure 
and return False, or 

(2) we have produced Vp. If it contains the identity, we return True, otherwise False. 

Notice that multiplying all elements of the set Vp by the inverse of any one of them 
produces the automorphism group of the block. 

3. Finding candidate Gram matrices 

We perform a backtracking search, based on the methods of Moyssiadis and Kou- 
nias |MKj and Chadjipantelis, Kounias and Moyssiadis jCKMj . to determine the set 
■M-n (d'^ain)- Starting from the 1x1 matrix, (n), candidate matrices are built up by 
symmetrically appending one row and column at a time, until size n is reached or no con- 
tinuation is possible. At this stage, the program returns to the most recent sub-matrix 
from which there is a possible continuation that has not yet been tried, and resumes the 
search from there. Proceeding exhaustively in this manner, until the entire search tree has 
been explored, the program terminates with a complete list of candidate Gram matrices. 

Naturally, the search tree is vast, and various methods must be used to prune it. The 
value d^in bounds the space we need to explore. In practice, a theorem of Moyssiadis and 
Kounias, described below, provides a criterion for removing branches of the search tree 
at an early stage. However, the size of the search space grows rapidly as dmin is lowered, 
so there is a practical limit to how low the threshold can be set. 
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A second method to reduce the size of the search is to impose lexicographic maximality 
of the matrices. This means that all but one of the branches containing a set of equivalent 
matrices are removed. There are two aspects to this pruning. The first is that only 
lexicographically ordered matrices whose blocks are in descending lexicographic order 
along the diagonal are generated. Specifically, there is an active block, namely the block 
containing the index of the most recently added row and column. A new block begins in 
index j whenever Pjj_i(M) is a minimal vector. The row and column vector appended to 
the current sub-matrix must have non-minimal elements only within the active block, and 
must be lexicographically less than the preceding row and column. (The first criterion is 
actually redundant.) Furthermore, a comparison is made of the active block (considered 
as a sub-matrix) and the previous block, and the continuation is disallowed if it would 
cause the blocks to be incorrectly ordered. 

However, as discussed in Section IT^ this alone does not guarantee that the matrix ob- 
tained is lexicographically maximal. To enforce maximality, the active block is subjected 
to the test IsLexMax, described in Section [2.31 and the continuation is disallowed if the 
test fails. This test can be extremely time-consuming, although it improved the efficiency 
of the algorithm for the cases considered here. For higher orders, it may be better to omit 
it, or, better, to carry it out more efficiently, perhaps using the methods of |McKj . 



3.1. The theorem of Moyssiadis and Kounias. 

Theorem 3.1. Let M = (^t be a symmetric, positive definite matrix of order 

m with elements taken from the set $ whose members are greater or equal in magnitude 
to some positive number c. Here and A are square matrices and the order of is 
r < m < n. The diagonal elements of A are equal to n. The columns of the r x {m — r) 
matrix B are taken from some set F,, C 



Define d* and 7* by 



d* 



D 



1 



r 
*T 



r 

c 



max 

7err 



Dr 

1' 



and define 



Then 



Ur{d) :- 



[n 



^ym-r-i ^(^^ _ + _ ^) max(0, d)] . 



detM < Ur{d*) 



(3.1) 



(3.2) 



(3.3) 



For the convenience of the reader, we will reproduce the proof of Moyssiadis and Kounias 
below. First we point out two slight modifications of their original statement of the 
theorem. 

• The columns of B and the vector 7 are taken to be members of F,. C $^ rather 
than of itself. This is needed in a practical implementation of the algorithm 
where lexicographic ordering is enforced. 
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• We do not require d* > but replace d by max(0, d) in the bound ()3.2|1 . This is 
necessary in the case n = 3 mod 4 since, in contrast to the case n = 1 mod 4, d* 
is often negative. 

Proof. By subdivision of the matrices A and B we rewrite M as 

'Dr 5i 7^ 
Bj A, 5 

where 7 G and (5 is a column vector of length m — r — 1 . Using the expansion by minors 
on the last column of M we obtain 

n R Dr Bi 
detM = {n-c) / + Bj A^ 6 



1 



6' 



Note that the first term is an (m — 1) x (m — 1) determinant that satisfies the hypotheses 
of the theorem. Denote the array in the second term by M. Either det M < or M is 
positive definite in which case we have 



det M = c 



Bl 



< C \Dr 



1 T I 

777 



The first factor is bounded above by d*, and the second, being the determinant of a 
positive definite matrix whose diagonal entries are less than or equal to tt, — c, is bounded 
above by (n — c)"^~^~'. Hence 



det M <{n 



\m—r—l 



max(0, d* 



By induction, using this bound and p.ip we obtain ()3.3|) . 



□ 



Corollary 3.2. Let M he an m x m symmetric, positive definite matrix whose diagonal 
elements are less than or equal to some positive number n and whose off-diagonal elements 
are greater than or equal to some positive number c in magnitude. Then 



det M < (n - c)"" + mc{n - c 



,m— 1 



Proof. Since M is positive definite, its determinant can only increase if we replace all 
diagonal elements with the upper bound, n. Now apply Theorem (|3.1|) with r = to this 
revised matrix. □ 

In the n = 1 mod 4 case our search algorithm will make use of Theorem ()3.1|) directly 
(with c = 1 and m = n). In the n = 3 mod 4 case we strengthen the theorem as follows: 



Corollary 3.3. Let M 



Dr B 

B^ A 



n = 3 mod 4 and furthermore let M E Ain,n- Then 

det M <{n- 1)"-" det Dr + [{n - l)"-" - {n - 3)"-'' -{n-r){n- 3)"""-^ 



satisfy the conditions of Theorem with m 

] max(0, d* 
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Proof. We apply Corollary 13 .21 to the matrix A[ := Ai — ^66^ in the proof of Theorem 13 .11 
to obtain a bound on detM which is sharper than ()3.1|) . Because all elements of both 
Ai and 6 are congruent to 3 modulo 4, it follows that the diagonal elements of A[ are 
bounded above hj n — 1 and the magnitude of the off-diagonal elements is bounded below 
by 2. Therefore 

det A[ <{n- 3^-'-^ + 2{m - r - l){n - 3^-'-^ 

and det M is bounded above by the product of this bound and max(0, d*). Applying the 
induction step, we find that the coefficient of max(0, d*) in the bound on det M becomes 
the sum 

n—1 — 1 

J2{n- - 3)^- + 2j{n - Sy^']. 

i=o 

Evaluating this sum gives the stated result. □ 



3.2. The algorithm. Let the set of allowed matrix elements be 

{l,-3,...,2-n} ifn = lmod4 
{—1, 3, ... ,2 — n} if n = 3 mod 4 



r : 


= 1 


Ml : 


= {n) 


Fo: 


= (()) 


ei : 


= n-2 


(Bei 


ginning 



and denote the determinant threshold (for {—1, l}-matrices) as d^nin- The threshold for 
candidate Gram matrices is then (i^in- 

We carry out the following recursive procedure. 

(1) Initialize variables. 

order of current sub-matrix 
the initial sub-matrix 

the list of vectors (Here () is the null vector.) 
the magnitude of the maximum allowed matrix element 

ihe recursive step) Increment r. Build the list of vectors, Fj—i, by 
appending to every vector in the list Fr-2 each e G $ satisfying |e| < e^ax- The 
construction is carried out in such a way that -Fr-i "will always be in ascending 
lexicographic order (according to the ordering <a). Initialize Fr-i to the null 
list. Build the list Tr by appending to every vector in Fr-2 all possible 2- vectors, 
(61,62) G satisfying |6i|, \e2\ < e^-i. 
(3) For each vector / in Fr-i construct the matrix 

li r = n then 

(a) If det Mr > and is a perfect square, then run the procedure IsLexMax 
on the final block of Mr. If it returns True then print out the matrix. 



THE MAXIMAL {-1, 1}-DETERMINANT OF ORDER 15 



11 



(b) Return to the calling procedure. 
If r < n go to step (4). 
(4) For each 7 G evaluate 

7 1 

until a (i is found such that the bound in Theorem 13.11 or Corollary 13.31 equals or 
exceeds d^^^, or until the list F,. is exhausted. If the list was exhausted without 
finding such a d, then is discarded, and we continue in step (3) with the 
next vector /. Otherwise, add / to the list Fj—i. Then test the active block for 
lexicographic maximality using IsLexMax. If it passes, recursively begin again at 
step (2), otherwise continue in step (3) with the next vector /. 
Notes: 

(a) If a new block started with index r, then index r + 1 is very important. It must be 
the case that all components of the vector f E that forms the initial segment 
of row and column r + 1 are minimal except possibly the r*^. If this component is 
minimal as well, then the entire remainder of the matrix M„ will be filled in with 
minimal elements, and we can proceed directly to the check of the perfect-square 
condition as in step (3). If component r is non-minimal, we set equal to its 
absolute value. 

(b) No attempt has been made to improve the efficiency of IsLexMax by keeping track 
of the automorphisms of the blocks of Dr as they are built. This may be necessary 
to extend to method to higher orders in the n = 3 mod 4 case. 

(c) No attempt has been made to improve the efficiency of step (4) by using facts about 
the form of the optimal form of 7. Moyssiadis and Kounias |MKj . for example, 
in their analysis of n = 17 were able to predict when a component would equal 
1 which avoided the need to check the entire list F^. This condition occurs much 
less frequently in the case n = 3 mod 4. 

4. Order by order list of candidate Gram matrices 

We state known results for all odd orders up to n = 21 and for selected higher orders. 
First we introduce a few standard types of matrices that arise repeatedly. 

Notation. The matrix Jp g is the p x q matrix with all elements equal to 1. We define 
Jp . Jp,p- 

Barba's bound is attained by the determinant of a { — 1, l}-matrix, R, if and only if 
Mr{R) = Mc{R) = {n — l)In + Jn- This occurs in orders 1, 5, 13, 25, and 41, among 
others, and we will have nothing more to say about these orders. No other candidate 
Gram matrix can have a determinant as large as Barba's bound, but when n = 1 mod 4 
it is not unusual to find candidate Gram matrices of high determinant which differ only 
slightly from the the matrix {n — l)In + Jn- Typically the off-diagonal entries differ from 
1 in only one or two rows and columns. 
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Notation. If f is a vector of lengtli n — 1 tlien S{v) denotes the matrix 



V 



{n - l)/„_l + Jn-l^ 

Notation. If v and w are vectors of lengtli n — 2 then D{a; v; w) denotes the matrix 




w 



l)/„_2 + J, 



n-2 , 



For n = 3 mod 4 define 
Definition. An Ehlich block matrix is a matrix of the form 



B{vi,V2, ...,Vk):={n- 3)/„ - + 4 



/J VI 









V2 



V 



\ 



with Y^Vj = n, or in other words, a matrix with diagonal elements equal to n, square 
blocks along the diagonal whose off-diagonal elements are equal to 3, and all other off- 
diagonal elements —1. 

In low orders, n = 3 mod 4, most large-determinant candidate Gram matrices are 
closely related to Ehlich block matrices. When an element a of a vector is repeated k 
times consecutively, it will be convenient to denote this sequence by a^. 

We describe below the candidate Gram matrices found by our procedure or, in a few 
cases which we note explicitly, matrices which were found by other means. In higher 
orders, where the the number of matrices tends to be large, we sometimes omit the 
explicit listing of matrices. Complete lists of matrices are available from the author upon 
request. 



4.1. n = 3. The largest determinant is 1 x 2^. Our program finds the unique candidate 
Gram matrix, -8(13), corresponding to this value, and no other candidates of equal or 
larger determinant. 

4.2. n = 7. The maximal determinant, found by Williamson |Wij . is 9x2^. It corresponds 
to the matrix 5(23, 1), which is found by our program. No other candidates of equal or 
larger determinant are found. 

4.3. n = 9. The maximal determinant, found by Ehlich and Zeller |EZj . is 56 x 2^. It 
corresponds to the candidate Gram matrix S{5, 17). No other candidates with equal or 
larger determinant are found. 
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4.4. n = 11. This case was treated by Ehlich, who found that the maximal determinant 
is 320 X 2^^, but his analysis was never published. His results were reported in the paper 
of Galil and Kiefer |GKj . There are seven candidate Gram matrices. Four of these have 
determinant (324 x 2^°)^, above the value corresponding to the maximal determinant. 
They are B{3, 2, Ig), 5(4, 17), B{5, U) and 



/II 3 3 
3 11 -1 
3 -1 11 



-J- 



3,8 



v 



8,3 



5(4, U 



J 



The failure of these four matrices to produce a maximal determinant is discussed in the 
next section. 

The other three candidates all correspond to the maximal determinant. They are 
5(5,23), 



/II 3 3 - - 

3 11 - 3 - 
3-11-3 

- 3 - 11 - 

- - 3 - 11 



-J- 



5,6 



V 



6,5 



and 



/II 3 3 3 - 

3 11 3 - - 

3 3 11 - - 

3 - - 11 3 

- - - 3 11 



-J' 



5,6 



6,5 



^(23) 



where "— " stands for —1. All candidates have unique characteristic equations. 



4.5. n = 15. All four candidate Gram matrices are Ehlich block matrices. The matrices 
5(43, 3) and 5(6, 3, 23) have determinant (105 ■ 3^ ■ 2^^^)^, corresponding to the conjectured 
maximal value, while 5(34,2,1) and 5(35) have determinant (108 ■ 3^ ■ 2^^)^. All of 
the matrices have different characteristic equations. Our algorithm, implemented as a 
Mathematica program, was able to generate the four candidates in about seven hours on 
a personal computer. This is a far larger running time than for any other order for which 
md(?T,) had been established, all of which require anywhere from a few milliseconds to a 
few minutes to produce the complete list of candidates. 



4.6. n = 17. The maximal determinant for n = 17 was conjectured by Schmidt jScj to be 
5 X 4^ X 2^^. The conjecture was independently formulated and proved by Moyssiadis and 
Kounias |MKj The candidate Gram matrices are S'(— 3, —3, 114) with determinant (22 x 
4^ X 2^y, S{5, -32, 113) with determinant (21 x 4^ x 2^^, 5(-3; -33, 112; I3, -3, In) with 
determinant (83 x 4*^ x 2^^f, 5(1; 5, -3, 113; 1, -33, In) with determinant (81 x 4*^ x 2^*^)^, 
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the matrix 





J2 


J2 


J2 


J2 




\ 


J2 


B 


J2 


J2 


J2 






J2 


J2 


B 


J2 


J2 






J2 


J2 


J2 


B 


J2 


Jl0,7 




J2 


J2 


J2 


J2 


B 










<^7,10 






I6J7 + J7 














/ 



where A 



17 5 
5 17 



B 



17 -3 
-3 17 



with determinant (5175 x 4 x 2 ) , and the six matrices 





= D(1;1, 


— Se, Is; 


—3, 114), 






= 5'(— 38, 


Is), 








= S'(— 3i6 


), 








= D(-3; 


—32, 113; 


—32, lis); 






= 5(52,- 


32, 112), 










/17 


-3 -3 


-3 1 


1 






-3 


17 -3 


1 1 


1 






-3 


-3 17 


1 1 


1 






-3 


1 1 


17 -3 


-3 




1 


1 1 


-3 17 


1 






1 


1 1 


-3 1 


17 



V 



11,6 



<^6,11 



16/n + J 



\ 



11 



each with determinant (5 x 4^ x 2^®)^. 

(17) 

The maximal determinant corresponds to the matrix M3 For some reason, the two 
matrices m'^'^^ and m'^'^^ were omitted from the hst in |MKj . The characteristic equations 
of the above matrices are distinct, except for those of and Mg '. 



4.7. n = 19. Our program fails to run to completion in a reasonable amount of time, or to 
turn up any candidate matrices, for any threshold below Ehlich's bound ()1.4|) . Two Gram 
matrices, corresponding to the best known determinant value, 833 x 4^ x 2^*, are found 
by direct computation from the conjectured maximal { — 1, 1} matrices of Smith |Smj and 
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Cohn |(]o2j . They are 



and 



/19 



19 

3 

3 



1 n 
19 
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o 


o 
O 


o 
O 










3 


19 


3 


3 
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— 




Q 
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Q 
O 


ly 


Q 
O 










Q 
O 


O 
O 


o 
o 


1 Q 

ly 


o 
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o 
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o 
O 


T 








o 
o 


1 Q 

ly 


o 
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o 
O 










o 
o 


o 
o 


1 Q 

ly 


Q 
O 










o 
O 


o 
o 


o 
o 


1 Q 

ly 


















5(34) 




3 






3 






3 


3 


3 














19 


3 














3 


19 


















19 


3 


3 












3 


19 


3 









10,9 



3 19 



19 

3 

3 



3 

19 

3 



3 
3 

19 



V 



9,10 



5(33 



Extensive searching by various random methods has failed to improve upon this deter- 
minant value. A simple exhaustive search for Ehlich block matrices with larger perfect 
square determinant than that of the above matrices yields the result that no such matrices 
exist. 

4.8. n = 21. The largest determinant for n = 21 is 29 x 5^ x 2^° which was found by 
Chadjipantelis, Kounias, and Moyssiadis |CKMj and proved by them to be maximal. The 
candidate Gram matrices are S'(— 35, 115) with determinant (6 x 5^° x 2^°)^, the matrix 

\ 





J2 


J2 


J2 


J2 


J2 


B 


J2 


J2 


J2 


J2 


J2 


B 


J2 


J2 


J2 


J2 


J2 


B 


J2 


J2 


J2 


J2 


J2 


B 



J^ 



10,11 



J 



11,10 



16/ 



11 



J 



11 



where A 



21 

5 



5 

21 



B 



21 

-3 



-3 
21 
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with determinant (18432 x 5^ x 2^°)^ and the three matrices Mf ' = 5(5,-35,114), 
Mf^^ = 5(5^1l6), and ikff^^ = 5(-7, -32, 117) all with determinant (29 x 5^ x 2'^°f. 

(21) (21) 

Ml ' and ' have the same characteristic equation; all others are distinct. The matrix 
corresponding to the maximal determinant is ■ For some reason the matrices M{ 
and Mf were omitted from the list in |CKMj . 

4.9. n = 29, 33, and 37. In orders 29 and 33 the largest known determinants were found 
by Bruce Solomon using the gradient ascent algorithm described in |OSDSj . In order 37 
the largest known determinant was constructed in |USj . 

Unfortunately, in all three cases these lower bounds are too low to use as thresholds in 
the search for candidate Gram matrices. Both time and space requirements appear to be 
growing exponentially as the threshold is lowered. The best we have been able to do is 
narrow the range within which md(n) must lie, by using somewhat higher thresholds. The 
thresholds used, and the number of candidates found for each determinant value equal to 
or exceeding the given threshold, are shown in Table ^ 

5. Decomposability of candidate Gram matrices 

5.1. General considerations. Given a candidate Gram matrix, M, we would like to 
determine whether it admits a decomposition M = RE7 where R is a { — 1, l}-matrix. 
This task appears to be rather difficult when only the matrix M is given. In the present 
context however, more information is available, namely the complete list of candidate 
Gram matrices with determinant equal to a specified perfect square value. In many cases 
one can demonstrate the nonexistence of decompositions using this knowledge. 

Let Mr be a candidate Gram matrix from our list, and suppose that = RR^ for 
some { — 1, l}-matrix R. Then, by Lemma 12.21 there is a matrix, = R^R, which has 
the same determinant and characteristic polynomial as Mj. Hence by permuting columns 
of R we can convert Mc to an equivalent matrix which also appears on our list. A first 
step in our procedure is therefore to enumerate the pairs {Mr,Mc) of candidate Gram 
matrices taken from our last which have the same characteristic polynomial. (Naturally, 
for every M on our list, (M, M) will be such a pair, but in general there may be other 
pairs too.) 

The pair of formulas 

Mr = RR^, = R^R (5.1) 

implies a number of constraints on the possible rows and columns of R, 

M^ = RM^R^ (5.2) 
M^ = R^MrR (5.3) 
Mr-R = RMc, (5.4) 

which are verified by substituting for Mr and Mc using ()5.1|) . Let xj denote the j^^ row 
of R, expressed as a row matrix, and let i/j denote the j^^ column of R, expressed as a 
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n = 29 


n = 33 


n = 37 


det. 


num. 


det. 


num. 


det. 


num. 


51 X 7 ^ '-^ 


2 


495 X 8^^ 


1 


680 X 9^^ 


1 


355 X 


1 


252315 X 8^^ 


1 


75 X 9^'' 


1 1 


352 X 7^2 


1 


3929 X 8^^ 


1 


672 X 9^^ 


1 


50 X 7^^ 


1 2 


490 X 8^'* 


1 1 


485070 X 9^^ 


1 


2448 X 7^^ 


1 


3919 X 8^^ 


1 


665 X 9^^ 


1 


119808 X 7^ 


1 


489 X 8^'* 


1 


53760 X 9^"^ 


1 


349 X 7^^ 


1 


3911 X 8^^ 


1 


4352000 X 9^^ 


1 


348 X 7^2 


1 


250047 X 8^^ 


1 


663 X 9^*^ 


1 1 


2432 X 7^^ 


1 


1023942465 x 8^ 


1 


483000 X 9^^ 


1 


2430 X 7^^ 


1 


3906 X 8^^ 


3 


4345600 X 9^^ 


1 


347 X 7^2 


2 


61 X 8^^ 


1-? 


5952 X 9^^ 


1 


118656 X 7^ 


1 


31203 X 8^^ 


1 


661 X 9^^ 


2 


2416 X 7^^ 


1 


3898 X 8^^ 


2 


5946 X 9^^ 


1 


345 X 7^2 


13 4 


3897 X 8^^ 


1 1 


53504 X 9^"^ 


1 


2413 X 7" 


1 


3895 X 8^3 


1 


5942 X 9^5 


2 


118128 X 7^ 


1 


31131 X 8^^ 


1 


660 X 9^*^ 


1-? 


344 X 7^2 


Is 


3889 X 8^=^ 


1,2 


659 X 9^*^ 


1 


2403 X 7" 


1 


31108 X 8^2 


1 






49 X 7^3 


l4,2 


486 X 8^^ 


15,3 






2400 X 7" 


l6,2 


31098 X 8^2 


1 






117504 X 7^ 


13 


3887 X 8^3 


1 






2397 X 7" 


3 


248724 X 8^1 


1 






5750784 X 7^ 


1 


3885 X 8^3 


1,2 






342 X 7^2 


13,2,3,7 


31050 X 8^2 
3881 X 8^3 


4 

1,1 







Table 1. Exhaustive enumeration of candidate Gram matrices of orders 
29, 33, and 37 with thresholds 342 x 7^^ 485 x 8^^, and 659 x 9^^ Square 
roots of the determinants of the candidates, with the factor 2^"^ omitted, 
are given. In number of matrices column, multiple entries refer to multiple 
characteristic equation classes. 



column matrix. The rows of Mr are denoted by mj, while those of Mc are denoted by 
m'J. Using the notation Xji for the i^^ element of the column matrix Xj, we have Xji = yij. 
The diagonal elements of ()5.2|) give a constraint on the rows of i?, 

mjmj = xjMcXj (5.5) 

while those of ()5.3p give a constraint on the columns, 

m'Jm'j = yjMryj. (5.6) 
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Finding solutions to these equations is facilitated by the large degree of symmetry that 
most candidate Gram matrices with large determinant have. A permutation of rows, 
represented by a permutation matrix P, is an automorphism of M if PMP^ = M. 
(Multiplying on the right by P"^ applies the same permutation to columns, preserving 
the symmetry of M.) The set of automorphisms forms a group under composition, the 
automorphism group of M. 

A special type of permutation, the swap, interchanges a pair of rows. The subset of the 
automorphism group generated by automorphisms which are swaps forms a subgroup of 
the full automorphism group. The permutation that interchanges rows i and j is written 
{ij). If (ij) and (ik) are automorphisms of M, then (jk) is as well. Hence we may order 
the rows of M so that the rows j such that (ij) is an automorphism, plus row i, form a 
contiguous block. The lexicographic ordering we imposed earlier already implies such a 
structure. It is then easy to see that the automorphism subgroup consisting of swaps is a 
direct sum of symmetric groups acting on these blocks. 

When such structure is present, the right-hand side of ()5.5|) will depend only on certain 
sums of elements of Xj, rather than on the elements individually. This greatly reduces the 
amount of work that needs to be done, as should become clear in the examples that follow. 
Note that we only generate solutions which are compatible with parity normalization (see 
Section 12)), that is, we require n + J2i ^ji = mod 4. Similar considerations apply to the 
equation for columns (j5.6j) . 

Once the row and column solutions have been found, we use the off-diagonal elements 
of ()5.2j) to determine which combinations of row (or column) solutions are consistent. The 
typical off-diagonal constraint takes the form 

mjnik = xjMcXk, (5.7) 

and involves the sums of elements of Xj mentioned above, as above, as well as the inner 
products of the sub-vectors of rows j and k corresponding to these sums. The analysis is 
most straightforward when the weights of these inner products are all equal. Then only 
the inner product of the full row, xjxk = rrijk, enters rather than the individual inner 
products. We will find that this occurs in all the examples studied below. 

There is also a compatibility condition between row and column solutions, derived 
from ()5.4|) . which is sometimes useful, namely 

mjyk = xjm'^. (5.8) 

In general the constraint will depend on the element i/kj = xjk and on the sums of elements 
referred to previously. In cases where the off-diagonal elements of the diagonal block of Mr 
containing rrijj and of the diagonal block of Mc containing m'^^ are equal, the dependence 
on ykj = Xjk is eliminated. In the examples explicitly worked out below, this constraint 
was not needed. 

The above constraints, sometimes combined with other elementary consequences of 
equation ()5.1|) . can often be used to prove the non-decomposability of a given pair of 
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candidate Gram matrices. We will illustrate the method on the candidate Gram matrices 
for order 15, and also on a certain pair of candidate Gram matrices of order 17. 

5.2. Order 15. The four candidate Gram matrices listed in Section E] correspond to 
different characteristic polynomials. Thus we must investigate four pairs 

(5(43, 3), 5(43, 3)), (5(6, 3, 23), 5(6, 3, 23)), 

(5(35), 5(35)), and (5(34, 2, 1), 5(34, 2, 1)). 

In each case we will write M = Mr = M^. Because the row and column solution sets are 
identical, some simplifications occur. 

5.2.1. (5(43, 3), 5(43, 3)). This Gram matrix corresponds to the conjectured maximal 
determinant. We will see that the decomposability analysis completely determines the 
form of the determinant. 

We may write M = 121 + 5 where 5 is the rank-4 matrix 

^ 3J4 — 1/4 — 1/4 — 

— J4 3J4 — J4 —^4,3 

— J4 — J4 3J4 —^4,3 
\ — ^3,4 —JsA —J3,4: 3J3 J 



B :-- 



Ifxisaj — 1,1} row or column vector, define 

4 8 12 15 

a := Xj b := Xj c := Xj d := Xj. 

j=l j=5 j=9 j=13 

It follows that a, b, and c are even and that d is odd, and that —4 < a,b,c < 4 and 
—3<d<3. Equations ()5.5p and ()5.6|) both reduce to 

83 = 4(a2 + b'^ + + d^) - (a + b + c + df 

for rows and columns 1 through 12. It is a simple matter to find the parity normalized 
solutions to this equation, subject to the above constraints, using an exhaustive search. 
The results are 

(a, b, c, d) = ([-4, -4, 0], 1), ([-4, -2, 2], -3), ([-2, -2, 2], 3), or ([-2, 0, 4], -1) 

where the notation [a, b, c] indicates that any permutation of a, b, and c is allowed. For 
rows and columns 13 through 15, the equation is 

75 = 4(a2 + b'^ + 0"^ + d^) - {a + b + c + df, 

and the solutions are 

(a, 6, c, d) = (-4, -4, -4, 1), ([-4, 0, 2], -1), (-2, -2, -2, 3), or ([-2, 2, 4], 1). 
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Not all of these solutions for, say, rows are compatible with each other. Imposing the 
equation (j5.7|) for all inequivalent choices of j and k reduces the solution sets to 

{(-4, -2, 2, -3), (-2, 0, 4, -1)} for rows 1-4, 

{(2, -4, -2, -3), (4, -2, 0, -1)} for rows 5-8, 

{(-2, 2, -4, -3), (0, 4, -2, -1)} for rows 9-12, 

{(-4, -4, -4, 1), (-2, -2, -2, 3)} for rows 13-15. 

Note that the first three of these sets may be permuted. Also the first three components 
of all solution vectors may be simultaneously permuted. Since such permutations give 
equivalent solutions, we disregard them. The column solutions must be also be compatible 
with each other, and also with the row solutions, according to equation ()5.8|) . The latter 
fixes a unique permutation of the column solutions, namely 

{(-4, 2, -2, -3), (-2, 4, 0, -1)} for columns 1-4, 

{(-2, -4, 2, -3), (0, -2, 4, -1)} for columns 5-8, 

{(2, -2, -4, -3), (4, 0, -2, -1)} for columns 9-12, 

{(-4, -4, -4, 1), (-2, -2, -2, 3)} for columns 13-15. 

For each block of rows (or columns) two different row (or column) types can occur. By 
using (j5.1|) we can determine precisely how many times each type occurs in the block. 
Consider rows 13 through 15. If any row has type (— 4, — 4, — 4, 1), that row will have 
inner product 11 or 15 with any other row of the same type. Hence a second row of type 
(—4, —4, —4, 1) is not allowed as the inner product must be 3. On the other hand, a type 
(—4, —4, —4, 1) row will have inner product 7 with any row of type (—2, —2, —2, 3), which 
is also disallowed. Therefore all three of the rows 13-15 must be of type (—2, —2, —2, 3). 
Now the matrix formed by the last three rows can be regarded as composed of three 3x4 
blocks followed by a 3 x 3 block consisting entirely of Is. The 3x4 blocks must have row 
sums —2 and column sums either —3 or —1. The row-sum condition implies that there is 
a single 1 in each row, while the column-sum condition implies that there can be no more 
than a single 1 in each column. Thus there is one column with column sum —3 and three 
with column sum —1 in each block. We may choose to place the —3 type column first in 
each block, and likewise for rows. Then the structure of the matrix is given by 

;-4, -2,2,-3) row 1 

-2,0,4,-1) rows2, 3, and4 

2, -4, -2, -3) row 5 

4,-2,0,-1) rows6, 7, and8 (5.9) 

-2,2,-4,-3) row 9 

0, 4, -2, -1) rows 10, 11, and 12 

-2, -2, -2, 3) rows 13, 14, and 15, 
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and 



-4,2,-2,-3) 

-2,4,0,-1) 

-2,-4,2,-3) 

0,-2,4,-1) 

2,-2,-4,-3) 

4,0,-2,-1) 

-2,-2,-2,3) 



column 1 

columns 2, 3, and 4 
column 5 

columns 6, 7, and 8 
column 9 

columns 10, 11, and 12 
columns 13, 14, and 15. 



(5.10) 



5.2.2. Uniqueness of the decomposition. The matrix found by Smith and Cohn exhibits 
the structure derived in the previous section. Here we will show that the matrix emerges 
quite simply from the structure and that there is essentially only one such matrix. 
The matrix must have the structure 



R 



Ml 
C 



A2 

C 
E2 



c 
B2 
A3 

E3 



DA 

D2 



(5.11) 



where Aj, Bj, and C are 4x4 sub-matrices; Dj are 4x3 sub-matrices; and Ej are 3x4 
sub-matrices. The matrices Aj have row and column sums (—4, —2, —2, —2); Bj have 
row and column sums (—2,0,0,0); and C has row and column sums (2,4,4,4) (which 
determines C completely and is why we need not distinguish different Cj). The matrices 
Dj have row sums (—3, —1, —1, —1) and column sums (—2, —2, —2), while for the Ej the 
row sums are (—2, —2, —2) and the column sums are (—3, —1, —1, —1). 
The row and column sums imply that 



/- 



Aj ~ A :-- 



A 



1 



/I 



B.r^B 



- 1 1 
1 - 1 
1 - 



D :-- 



1 



and 



1 



E, ~ E : - 



C 



(- 
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1\ 


1 


1 


1 


1 


1 


1 


1 


1 


V 


1 


1 


V 


1 - 


) 






- 1 









where "~" means equivalent up to permutation of the last three rows or columns. The 
structure ()5.9|) . ()5.10|) . ()5.1H) is preserved under permutation of rows 2-4, of rows 6-8, 
of rows 10-12, of rows 13-15, or of the corresponding columns. We use the freedom to 
permute columns to fix Ai = A, Bi = B, and Di = D in ()5.11|) . Then use the freedom to 
permute rows to fix A2 = A, B^ = B, and Ei = E. Finally, another column permutation 
fixes B2 = B. Thus rows 1-4 and columns 1-4 are completely constructed and have 
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the correct pairwise inner products. In order for these to have the right inner products 
with rows 5-8 and columns 5-8, we must have D2 = D and E2 = E. For rows 13-15 to 
have the desired inner products with rows 6-8 necessitates = E. Similarly columns 
13-15 have the correct inner products with columns 2-4 only if E^ = E. It then follows 
that A3 = A, and the matrix R such RR^ = R^R = M is completely determined up 
to permutations consistent with the structure of M, and negation of the entire matrix. 
Consequently there is essentially only one D-optimal design of order 15 (subject to the 
proof of non-decomposability of the other three candidate Gram matrices, which follows). 

5.2.3. (-8(6, 3, 23), -8(6, 3, 23)). The candidate matrix may be written M = 121 + B where 
-B is a rank-5 matrix. Row and column vectors are subdivided in to 5 sub-vectors whose 
element sums are denoted a, b, c, d, and e. The sum b is odd whereas the other four are 
even. The inequalities —6 < a < 6, —3 < 6 < 3, and —2 < c,d,e < 2 must hold. 
For rows or columns 1 through 6, equations ()5.5|) and ()5.6|) reduce to 

99 = 4(a^ + b'^ + c^ + d^ + e^)-{a + b + c + d + ef. 

Exhaustive search yields the solution set 

(a,6,c,rf,e)G{(-6,-l,0,0,0), (-4, 1, [0, 2, 2]), (-4,3,-2,-2,-2), 

(-2,-3,2,2,2), (2, -3, [-2, 2, 2]), (4,-3,0,0,0), 

(4,1, [-2, -2,0]), (6,-1, [0,2, 2]), (6,3,0,0,0)}. 

For rows or columns 7 through 9, equations ()5.5|) and ()5.6|) reduce to 
75 = 4(a^ + b'^ + + d^ + e^) - {a + b + c + d + ef 
which has the solution set 

(a, 6, c, rf, e) G {(-6, -3, [-2, 0, 0]), (-6, 1, -2, -2, -2), (-4, -1, [0, 0, 2]), 
(-2, 3, [-2, -2,0]), (0,3,-2,-2,-2), (2,-3,2,2,2), 
(4,1, [-2, 0,2])}. 

For rows or columns 10 through 15, equations ()5.5|) and ()5.6|) reduce to 

67 = 4(0^ + b'^ + + d^ + e^)-{a + b + c + d + ef 

which has the solution set 

(a, 6, c, rf, e) G {(-4, -1, [-2, -2, 2]), (-2, 1, [-2, 2, 2]), (0, -3, [0, 2, 2]), 

(2, -3, [0,0, 2]), (2,1, [-2, -2, 2]), (4, 3, [-2, 2, 2]). 

If for a given row 1-6 solution the compatibility condition (|5.7|) fails to hold for all row 
7-9 solutions, we drop that solution. Likewise if a given row 7-9 solution is incompatible 
with all row 1-6 solutions, we do not consider that solution any further. The compatibility 
condition in this case reduces to 

—33 = (4/5 — J^)xi 
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where Xj is a solution for row j taken from the appropriate hst above. Failure to find any 
such that this equation holds when xi is one of (—4,1, [0,2,2]), (—4, 3, —2, —2, —2), 
(—2,-3,2,2,2), or (4, 1, [— 2, — 2, 0]) eliminates those solutions. Similarly the solutions 
X7 = (-2, 3, [-2, -2, oj) and (4, 1, [-2, 0, 2]) are eliminated. 

Imposing compatibility of the remaining row 1-6 solutions with the row 10-15 solutions 
allows us to eliminate the solutions xio = (0,-3, [0,2,2]) and (2,1, [-2,-2,2]), and at 
the same time, the solutions xi = (4, —3, 0, 0, 0) and (6, 3, 0, 0, 0). As a consequence, xr = 
(2, —3, 2, 2, 2) is also eliminated, since the only xi solution with which it was compatible 
was (6,3,0,0,0). 

Compatibility of the row 7-9 solutions with the row 10-15 solutions then eliminates 
Xio = (-2,l,[-2,2,2]). 

Finally, one may check that the compatibility condition between, say, row 10 and row 
12, eliminates all the remaining row 10-15 solutions as possibilities. Therefore, there is 
no decomposition of (5(6, 3, 2, 2, 2), B{6, 3, 2, 2, 2)). 

5.2.4. (5(3, 3, 3, 3, 3), 5(3, 3, 3, 3, 3)). The equations (j5.5|) and (j5.6|) for the sums of row 
and column sub-vectors have the same form for all 15 rows and all 15 columns, namely 

75 = 4(a^ + 1"^ + + (f + e"^) - {a + b + c + d + ef , 

where a, 6, c, d, and e are odd integers in the range [—3,3]. The parity normal- 
ized solutions to this equation are {a,b,c,d,e) = ([-3,-3,1,1,1]) and {a,b,c,d,e) = 
([—3, —1, —1, —1, 3]). Since all solutions have at least one variable set to —3, we may use 
our freedom to permute columns to fix a = —3 in the row 1 solution. Likewise, we fix 
b = ±3 in row 1. Now if {a,b,c,d,e) = (-3,-3,1,1,1) in row 1, then the consistency 
equation (|5.7|) for rows requires that the same solution hold for rows 2 and 3 and that 
each of rows 4-15 be one two types: (a, b, c, d, e) = (—3, 3, —1, —1, —1) or (a, b, c, d, e) = 
(3,-3,-1,-1,-1). On the other hand, if {a,b,c,d,e) = (— 3, 3, — 1, — 1, — 1) in row 1, 
then the same solution must hold for rows 2 and 3, and the solution (a, b, c, d, e) = 
(—3, —3, 1, 1, 1) must obtain in rows 4-15. 

The combination of these two facts and permutation symmetry imposes the requirement 
on the five sets of rows, {1, 2, 3}, {4, 5, 6}, . . ., {13, 14, 15}, that within each set all three 
rows must be of the same type, but that outside the set, all rows must be of the other 
type. As no more than two row sets may be constructed under these constraints, whereas 
we require five, no decomposition exists. 

5.2.5. (5(3, 3, 3, 3, 2, 1), 5(3, 3, 3, 3, 2, 1)). The impossibility of decomposition follows by 
an argument that is similar to but more involved than that of the previous case. We need 
only look at the solution set for rows 1-12. Equation (j5.5j) becomes 



75 = 4(a2 + b^ + c'^ + d^ + e'^ + f)-ia + b + c + d + e + ff 
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where a, b, c, and d are odd integers in the range [—3, 3], e is an even integer in the range 
[—2, 2], and / G {—1,1}. The parity normahzed solutions are in the set 

{([-1,-1,3,3], 2,-1), ([-3,1, 1,3], 2,1), ([1,1,3,3], -2,-1), 
([-1,1, 3, 3], -2,1), ([-3, -3, 1,1], 0,1), ([-3, -1,-1, 3], 0,-1)}. 

We now use the compatibihty conditions among rows ()5.7j) to prove the impossibihty 
of a decomposition. We leave the arithmetical checking of the consistency of various pairs 
of solutions to the reader. A requirement of decomposability is that we be able to assign 
consistent solutions to the four sets of rows {1, 2, 3}, . . ., {10, 11, 12}. If the set {1, 2, 3} 
has a row of type (—1, —1, 3, 3, 2, —1) then rows 4-12 must be of type (1, 1, [—3, 3], 2, 1). 
But two rows, in different row sets, which both are of the latter type are not compatible 
with each other, so it is not possible to construct four sets when one of them contains a 
row of type ([—1, —1, 3, 3], 2, —1). 

Similarly, if the set {1,2,3} contains a row of type (—3,1,1,3,2,1), then rows 4-12 
must all be of type (1, 3, 3, 1, —2, —1). Again, two rows of this type in different row sets 
are not compatible, so (—3, 1, 1, 3, 2, 1) has been eliminated. 

If either of the two types (—1, 1, 3, 3, —2, 1) and (1, 1, 3, 3, —2, —1) appears in the set of 
rows {1, 2, 3}, then rows 4-12 are of a suitable form of type ([—3, —3, 1, 1], 0, 1) or type 
([—3, —1, —1, 3], 0, —1). This means that at least two of the three sets {4, 5, 6}, {7, 8, 9}, 
{10, 11, 12} must both contain a row of the first type, or both contain a row of the 
second type. However neither of the latter two types can appear in two different row sets. 
Hence the solutions ([—1, 1, 3, 3], —2, 1) and ([1, 1, 3, 3], —2, —1) are ruled out. But since, 
as we just argued, neither of the two remaining solution types, ([—3, —3, 1, 1],0, 1) and 
([—3, —1, —1, 3], 0, —1) can appear in two different row sets, constructing four row sets is 
impossible. Therefore there is no decomposition. 

5.3. Other orders. In orders 3, 7, and 9, decompositions of the the unique largest can- 
didate Gram matrices have been known for some time. 

5.3.1. n = 11. Ehlich proved the non-decomposability of the four candidates of largest 
determinant, and showed that the three candidates of the next largest determinant do 
decompose. We have applied our method to these seven matrices, and confirmed his 
results. As the method does not differ in any substantial way from that applied to the 
case n = 15, we omit the details. 

5.3.2. 77, = 17. In |MKj eight of the nine reported candidate Gram matrices were shown 
to be non-decomposable, thus establishing the maximality of the determinant associated 
with the ninth. Since our list contains two additional candidates, we must, in order 
to verify that result, prove them non-decomposable as well. Furthermore, one of the 
additional candidates has the same characteristic equation as one of the original nine, so 
this pair must also be checked. 
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The three decompositions to be tested are then 

{S{52, —32, 112), 5'(52, —32, 112)), 

{D{1- 1, -36, Ig; 3, Im), D{1; 1, -36, Ig; 3, l^)), 

(5(52,-32,ll2),5(-38,l8)). 



For the former two, follow the model of proof in |MKj . (We have checked these, as well as 
the original nine candidates.) As the last pair gives us a chance to illustrate the method 
when Mj. 7^ M^, we present the proof here. 

We let Mr = S{52, —32, 112) and Mc = 5'(— 38, Is)- We compute the sets of solutions to 
the quadratic equation ()5.5|) for row 1 and for rows 2-3, and then show that the sets are 
incompatible. The equation for row 1 is 

97 = + 6^ + - Gab + lac + 26c 

where a G {—1,1} and h and c are even integers in the range [—8, 8]. The parity normalized 
solutions are (a, 6, c) = (1, —2, —8) and (a, 6, c) = (—1, 6, 2). For rows 2 and 3, the equation 
is 

57 = + 6^ + - <oah + 2ac + 26c 

and the parity normalized solutions are (-1,4,-8), (-1,6,-2), (1,-6,8), (1,-4,-2), 
and (1,8,2). The compatibility constraint (|5.7|) 

mi ■ m2 — 16mi 2 = x\{Mc — 16/)x2 

which reduces to 

101 = aia2 + 6162 + C1C2 - 3(0162 + 6102) + ctiC2 + Cia2 + 61C2 + C162 

where x\ is represented by (ai,6i,ci) and by (a2,62,C2). Taking all pairs of solutions 
from the sets enumerated above, we find that none satisfy the constraint. Therefore no 
decomposition exists. 

5.3.3. n = 19. There are presently no known candidate Gram matrices. 



5.3.4. n = 21. As in the case n = 17, in order to verify the result of |CKMj we have three 
additional pairs of Gram matrices to prove non-decomposable, 

{S{5, —35, 114), S{5, —35, 114), 

(S'(-7, -32, 117), 5'(-7, -32, 117)), and 

(5(5,-35,ll4),5(-7,-32,ll7)). 

The proof parallels that of the n = 17 case, so we omit it here. 
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5.3.5. n = 29, 33, and 37. All of the candidate Gram matrices enumerated in Table C] 
have been shown not to decompose. As there are a great many matrices to check, an 
automated procedure was used for this purpose which combines some, but not all, aspects 
of the method illustrated in the case of n = 15 with a backtracking search. (In fact, the 
result for n = 15 was originally obtained by the use of this algorithm. As its performance 
is far from optimal, we will not describe it in detail here. We hope to present an improved 
version in a forthcoming publication.) 

Combining the above with the results of |USDSj and jUSj we are able to place new 
bounds on md(n) for these three orders. For = 29 we conclude that 320 x 7^^ < md(29) < 
342 X 7^2 or in terms of Barba's bound (HID, 0.865001 . . . < md(29)/5(29) < 0.92447. 
For n = 33 we have 441 x 8^^ < i^(33) < 485 x 8^^ or 0.854667. . . < md(33)/5(33) < 
0.939951. Finally, for n = 37 we find 72 x 9^^ < ^{37) < 659 x 9^^ or 0.936329. . . < 
md(37)/fi(37) < 0.952224. 



6. The range of the determinant function for n = 9 and n= 11 

Another application of these methods is the determination of the complete range of 
the determinant function for { — 1, l}-matrices. By negation of appropriate rows and 
columns, we can make the first row of a matrix consist entirely of Is and the first column, 
except for the (1, 1) entry consist entirely of —Is. Such matrices are called normalized. 
Adding the first row of a normalized matrix to every other row does not change the 
determinant, and produces a matrix whose first column has a single 1 in the (1,1) position 
and Os elsewhere. All entries in rows 2 through n, where n is the size of the matrix, 
are either or 2. Expansion of the determinant by minors on the first column then 
shows that the determinant is equal to that of a {0, 2}-matrix of size n — 1. Thus an 
n X n { — 1, 1 }-determinant is always divisible by 2"~^. Furthermore, we have shown the 
problem of maximizing the determinant of n x n { — 1, l}-matrices to be equivalent to that 
of maximizing {n — 1) x [n — 1) {0, l}-matrices. 

The problem of the range of the determinant function was studied by Craigen |Crj 
who asked for the complete list of integers, d, such that d is the determinant of some 
(n — 1) X (n — 1) {0, l}-matrix, or equivalently, 2"'~^d is the determinant of some n x n 
{ — 1, l}-matrix. We are focusing here on { — 1, l}-matrices, but it is convenient to omit 
the factor 2"~^ in the discussion, so we will be careful always to refer to normalized 
determinants (see Section H}. 

We ignore the sign of the determinant. For n < 7 the range of the normalized deter- 
minant function is [0,md(n)], where the notation [a,b] means all integer values between 
a and b. Craigen proved that in n = 8, a gap appears for the first time. In particular he 
showed that there are no normalized determinants in the range [28, 31] whereas the order 
8 Hadamard matrix has normalized determinant 32. 

Craigen also provided a table of normalized determinant values known to him at the 
time. For n = 8 and n = 9 we reproduce his lists in Table El (For n = 8 Craigen also lists 
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n = 8 


n = 9 


n = 10 


n = 11 


\0 181 


\0 401 


fO 1021 


\0 2681 


20 


42 


104 1051 


270 2761 


94 


[44 4^;! 






32 


48 


110 


[282, 286] 




56 


112 


288 






[116, 117] 


291 






120 


[294, 297] 






125 


304 






128 


312 






144 


315 








320 



Table 2. The range of the normahzed determinant function for orders 8, 
9, 10, and 11. 



19, which is almost surely a misprint.) For orders 10 and 11, Craigen gave partial lists. 
In the table, we give what we believe to be complete lists of normalized determinants. 

We have been able to exhibit a determinant for each of the values in the table. The 
method by which this was done is beyond the scope of this paper, and we intend to 
describe it elsewhere. What is of relevance to us here is that for the two odd orders, 
n = 9 and n = 11, we are able to prove that the lists are complete, that is, there are no 
normalized determinant values other than those given in the table. 

The method is as follows: we set the threshold in our candidate Gram matrix finding 
routine to the lowest gap in the range. For n = 9 this is 41 x 2^ and for = 11 it 
is 269 X 2^°. The program then produces a complete list of candidate Gram matrices 
corresponding to determinants equaling or exceeding the threshold. If none of these 
determinants correspond to any of the gaps, then we have proved the claim. On the other 
hand, if any of the candidate Gram matrices do correspond to gaps, we must then show 
that these matrices do not decompose. 

In the case n = 9 we find candidate Gram matrices corresponding to the following 
normalized determinants: 56 (1 matrix), 48 (4 matrices), 45 (1 matrix), 44 (2 matrices), 
and 42 (1 matrix). Since none of these correspond to gaps, completeness is proved. 

For n = 11 we find matrices corresponding to the determinants listed in Table El There 
are 196 matrices in total. All of the determinant values correspond to gaps except for 
three of them, namely 300, 306, and 324. The non-decomposability of the four matrices 
with determinant (324 x 2^°)^ was discussed previously. To establish the desired result, 
one need only show that the three matrices with determinant (300 x 2^°)^ and the matrix 
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det. 
mult. 


270 
47 


271 
1 


272 
37 


273 

8 


274 
2 


275 
8 


276 
5 


278 
1 


279 
7 


280 
10 


det. 
mult. 


282 
2 


283 
1 


284 
2 


285 
2 


286 
1 


288 
41 


291 
1 


294 
2 


295 
1 


296 
1 


det. 
mult. 


297 
1 


300 
3 


304 
2 


306 
1 


312 
1 


315 
1 


320 
3 


324 
4 







Table 3. Multiplicities of the candidate Gram matrices corresponding to 
the listed normalized determinants of order 11. The threshold used is 269. 



with determinant (306 x 2 ) do not decompose. These matrices are 







3 


3 








3 


11 




3 






3 




11 




3 —^5,6 






3 




11 


3 








3 


3 


11 










,5 


S(32) 
















3 


3 






\ 


3 


11 




3 






3 




11 




3 






3 




11 


3 








3 


3 


11 
















v 










J 



\ 



/II 3 3 - 

3 11 - 3 

3 - 11 3 

- 3 3 11 



-J. 



4,7 



-J: 



7,4 



i?(32,l) 



and 



/II 3 3 3 - 

3 11 - - 3 

3 - 11 - - 

3 - - 11 - 

- 3 - - 11 



-J' 



5,6 



-Ji 



\ 



6,5 



S(3,l3 



There are no unusual features in the proofs of non-decomposability, so we omit them. 

7. Conclusion and outlook 

We have proved that the maximal determinant of a 15 x 15 { — 1, l}-matrix is 25515 x 2^^ 
and confirmed the known maximal values for the other odd orders which do not attain the 
Barba bound, up to order 21, except for 19. We have also established new upper bounds 
on the maximal determinant in orders 29, 33, and 37. 

Convincing conjectures exist for the maximal determinants for orders n = 19 |Smj . 22, 
23 OSDSJ, and 37 |USj . We hope that the methods of this paper can be extended to 
handle at least some of these cases. 

Order 19 may be tractable using the current method with efficiency improvements in 
the computer code, and perhaps parallelization. For higher orders, n = 3 mod 4, there are 
excellent lower bounds on the maximal determinant, but these can possibly be improved. 
It is unlikely that our method can be applied to these cases without a new algorithm. 
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The smallest even order for which the maximal determinant is unknown is 22, which 
may well be tractable, although it is too early to tell until experiments are done. 

The biggest hope is for orders n = 1 mod 4, the lowest open case of which is n = 29. 
The major hurdle to be overcome here is the inefficiency of the program that searches for 
candidate Gram matrices. Efforts are underway to try to improve this. In general, the 
n = 1 mod 4 cases are the most amenable to our method, and the major determinant of 
success is how close the actual maximum lies to the Barba bound. Thus we are especially 
optimistic about the case n — 37. 
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